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^ Abstract 

The thermochemistry of the carbon clusters C n (n=2-10) has been revisited by means of W4 
theory and W3.21ite theory. Particularly the larger clusters exhibit very pronounced post-CCSD(T) 

OS 

correlation effects. Despite this, our best calculated total atomization energies agree surprisingly 

well with 1991 estimates obtained from scaled CCD(ST)/6-31G* data. Accurately reproducing 

the small singlet-triplet splitting in C2 requires inclusion of connected quintuple and sextuple 
g 

excitations. Post-CCSD(T) correlation effects in C4 stabilize the linear form. Linear/cyclic 

q equilibria in Cq, Cg, and C10 are not strongly affected by connected quadruples, but they are 

q affected by higher-order triples, which favor polyacetylenic rings but disfavor cumulenic ones. 

• 1— 1 

Near the CCSD(T) basis set limit, C10 does undergo bond angle alternation in the bottom-of- 

43 

the-well structure, although it is expected to be absent in the vibrationally averaged structure. 
The thermochemistry of these systems, and particularly the longer linear chains, is a particularly 
difficult test for density functional methods. Particularly for the smaller chains and the rings, 
qq double- hybrid functionals clearly outperform convential DFT functionals for these systems. Among 

compound thermochemistry schemes, G4 clearly outperforms the other members of the Gn family. 
Our best estimates for total atomization energies at K are: C2( 1 S+) 144.07, C2( 3 n u ) 142.39, 
Csi 1 ^) 315.83, C 4 ( 3 S 9 ) 429.16, C^A g ) 430.09, C 5 ( 1 S+) 596.64, C 6 ( 3 £ g ) 717.19, C^A'J 
729.68, C 7 ( 1 S+) 877.45, C 8 ( 3 £ 9 ) 1001.86, C 8 (M 9 ) 1014.97, C 9 ( 1 S+) 1159.21, Ci ( 3 X 9 ) 1288.22, 
S and Cio^Ai) 1355.54 kcal/mol. 
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I. INTRODUCTION 



Carbon clusters, C n , are on the one hand precursors for diamond films generated by CVD 
(chemical vapor deposition), and on the other hand the smaller congeners of fullerenes and 
carbon nanotubes. Their chemistry has been reviewed, among others, by Van Orden and 
Saykally[lJ, by Lifshitz[2j, and most recently in the Introduction to Ref.[3]. 

One of the earliest papers to offer a coherent picture of chemical bonding in the C n clusters 
for lower n was the seminal study of Raghavachari and Binkley (RB) [4J. While based on what 
by present-day standards are fairly primitive calculations (HF/6-31G* optimizations and 
frequency calculations accompanied by CCD(ST)/6-31G* single-point energy evaluations), 
RB established a number of key features of the chemistry of these systems. This work 
was expanded upon to some degree by Hutter, Liithi, and Diederich (HLD) [5] , who carried 
out DFT optimizations and frequency calculations through Cis, and by Martin and Taylor 
(MT1)[6], who carried out CCSD(T) coupled cluster optimizations and harmonic frequency 
calculations with fairly modest basis sets, reviewed the experimental information available 
at this point, and proposed a number of reassignments of infrared spectroscopic features. 
Guided by predictions in this latter study, cyclic Cq^\ and Cs[E] were experimentally 
discovered in graphite vapor trapped in a solid argon matrix. (The first experimental 
evidence for nonlinear small carbon clusters were probably Coulomb explosion experiments 
carried out at the Weizmann Institute [U].) 

Earlier, linear C4 and C5 were correctly assigned[10j following computational 
predictions [TT], and the characterization of linear C6 likewise involved an interplay 
between theoryp.2] and experiment |13j. Linear C7 and C9 were discovered [H] following 
computational predictions |12j. 

As discussed by RB and MT1, the odd- numbered clusters C2«+i have linear structures 
with closed-shell 1 S+ ground states, while for the even-numbered clusters, linear structures 
in 3 S ff ground states are energetically close to closed-shell singlet ring structures. Among 
the latter, the 64^+2 rings are cumulenic with D^ n +i)h symmetry and exhibit twin aromatic 
systems (one conventional, another in-plane): while the ring strain on the in-plane system 
in C6 is simply too great and leads to both alternating-angle distortion and the ring being 
less stable than the chain. C10 is the first carbon cluster for which the ring structure is more 
stable than the chain. As for the C^n clusters, while the diamond-shaped C4 ring is nearly 
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isoenergetic with the chain[T5l [16], the Cg ring is polyacetylenic and exhibits both bond 
length and bond angle alternation (C^ symmetry), and larger C 4n rings were predicted [TT] 
to generally have polyacetylenic C^ n )h structures. At some point for high 4n + 2, Peierls 
distortion will set in: this issue has been discussed and reviewed at length in Ref . [3] . 

Thermochemistry was addressed in a number of theoretical and experimental studies. RB 
applied an empirical scaling factor of 1.1 to their CCD(ST)/6-31G* atomization energies, 
which they derived from the ratio between their calculated values and Knudsen effusion 
measurements by Drowart et al.(DBDI)[18j for C2 through C5. Since the DBDI values for 
C4 and C5 were third-law extrapolations based on very crude estimates of the molecular 
constants, the scaled estimates of RB were biased upwards, and the scaling factor required 
downward revision to 1.082 [19]. Later, new measurements by Gingerich and coworkers|20] 
became available, which also covered C6 and C7. Subsequently, Martin and Taylor 
(MT2)[2I] proposed revised estimates based on CCSD(T)/cc-pVTZ and CCSD(T)/cc-pVQZ 
calculations with empirical basis set incompleteness corrections. 

Both computational chemistry in general, and computational thermochemistry in 
particular, have come a long way since this time. In particular, our group recently developed 
a post-CCSD(T) computational thermochemistry protocol known as W4 theory [22]. which 
exhibits an RMSD deviation from experiment (Active Thermochemical Tables, ATcT[23j, 
values) of just 0.08 kcal/mol, implying a 95% confidence interval of just 0.16 kcal/mol. The 
purpose of the present paper is to furnish the best thermochemistry for the lower C n clusters 
feasible with current technology. 

Carbon clusters are among the manifold research interests of Prof. Henry F. Schaefer 
III [21] EH] [26] ET] [28], who is being honored by the present issue. 

II. COMPUTATIONAL METHODS 

The self consistent field (SCF), ROCCSD and ROCCSD(T) calculations [29] were carried 
out using version 2006.1 of the Molpro program system [30J. All single-point post- 
CCSD(T) calculations were carried out using an OpenMP-parallel version of Mihaly Kallay's 
general coupled cluster code MRCCjM] interfaced to the Austin-Mainz-Budapest version of 
the ACES II program system[32j. The diagonal Born-Oppenheimer correction (DBOC) 
calculations were carried out using the Austin-Mainz-Budapest version of the ACES II 
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program system as well as using PSI3[33J. 



For the lower-cost methods W2.2 and W3.21ite we used reference geometries optimized 
at the B3LYP/cc-pVTZ level of theory, while for the more rigorous methods W3.2, W41ite, 
and W4 we used CCSD(T)/cc-pVQZ reference geometries. The optimized geometries are 
summarized in Table [H 

For the large scale SCF, CCSD, CCSD(T), and post-CCSD(T) single point calculations, 
we employed the cc-pVnZ basis set[34j. In core- valence correlation calculations, the 
core-valence weighted correlation consistent basis sets of Peterson and Dunning were 
employed[35j. Scalar relativistic calculations were carried out using the Pacific Northwest 
National Laboratory (PNNL) Douglas-Kroll-Hess relativistically contracted correlation basis 
sets 1 361. 



The ROHF-SCF contribution is extrapolated using the Karton-Martin[37] modification 
of Jensen's extrapolation formula. |38j All other extrapolations are carried out using the 
A + B/L a two-point extrapolation formula (where L is the highest angular momentum 
present in the basis set). The ROCCSD valence contribution is partitioned into singlet-pair 
energies, triplet-pair energies and T\ terms. [39] The singlet- and triplet-pair energies are 
extrapolated with as =3 and =5, respectively, while the T\ term (which exhibits very 
weak basis set dependence) is simply set equal to that in the largest basis set. All other 
extrapolations are carried out with a =3. [221 HO] 

The Wn family of methods Wl, W2.2, W3.21ite, W3.2, W41ite, and W4 used in the 
present study provide a sequence of converging computational thermochemistry protocols. 
A detailed description and rationalization of the Wn protocols is given elsewhere. [221 |4"T1 l4*2l 
H3J 01] For the purpose of the present paper we use the Wnh variants of the Wn methods, 
in which the diffuse functions are omitted from carbon and less electronegative elements[42j. 
(For the present all-carbon systems, this approximation is of no thermochemical consequence, 
but computer resource requirements are substantially reduced.) In W4h theory the SCF and 
valence CCSD contributions to the TAE are extrapolated from cc-pV5Z and cc-pV6Z basis 
sets, and the valence parenthetical triples (T) contribution from cc-pVQZ and cc-pV5Z 
basis sets. The higher-order connected triples, T3 — (T), valence correlation contribution is 
extrapolated from the cc-pVDZ and cc-pVTZ basis sets. As for the connected quadruple, T 4 , 
term — the (Q) and T 4 — (Q) corrections are calculated with the cc-pVTZ and cc-pVDZ basis 
sets, respectively, both scaled by 1.1. This formula offers a very reliable as well as fairly 




4 



cost-effective estimate of the basis set limit T4 contribution. j22j HO] The T5 contribution 
is calculated using the sp part of the cc-pVDZ basis set, denoted cc-pVDZ(no d). The 
CCSD(T) inner-shell contributions are extrapolated from cc-pwCVnZ basis sets (n =T, Q). 
Scalar relativistic contributions (second-order Douglas-Kroll-Hess approximation [IS]) are 
obtained from the difference between nonrelativistic CCSD(T)/cc-pVQZ and CCSD(T)/cc- 
pVQZ-DK calculations. Atomic and molecular first-order spin-orbit coupling terms are taken 
from the experimental fine structure. Finally, the DBOC is calculated at the ROHF/cc- 
pVTZ level of theory. 

When considering the Wn sequence, it may be helpful to keep the following points in mind: 
(a) Wl and W2.2 completely neglect post-CCSD(T) correlation effects; (b) the difference 
between Wl and W2.2 is (in the context of this paper) principally a matter of more reliable 
CCSD(T) basis set limits and the DBOC; (c) the only differences between W2.2 on the one 
hand and W3.21ite and W3.2 on the other hand are post-CCSD(T) corrections, the former 
with a cost-saving empirical approximation; (d) the only difference between W3.2 and W41ite 
is a more reliable basis set limit for the CCSD(T) part; (e) the only improvement in W4 
relative to W41ite is a more rigorous account for connected quadruple and higher excitations. 

All density functional results were obtained by means of a locally modified version of 
Gaussian 03[46j. 



III. RESULTS AND DISCUSSION 

Diagnostics for the importance of nondynamical correlation can be found in Table [TT} 
These include the T\ and T>\ diagnostics [4TI 148] . the largest T 2 amplitudes at the CCSD/cc- 
pVTZ level, the HOMO and LUMO natural orbital occupations at the same level, and 
percentages of the total atomization energy accounted for at the SCF level, by parenthetical 
triples (shown elsewhere [22] to be a more reliable predictor of the importance of post- 
CCSD(T) correlation effects than any other diagnostic), and by post-CCSD(T) correlation 
effects. 



Table |III| presents component breakdowns for the W3.21ite results, Table |IV| for the 
W4 results (C2 through C5 only), and Table [V] a comparison between theoretical and 
experimental estimates for C2 through Cio- 

For C2, a well-known 'problem molecule', [23 SSI ED] an ATcT value is available: [5T] our 
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W4 calculation agrees with it to within overlapping uncertainties. Still higher-level W4.3 and 
W4.4 calculated values are in basically perfect agreement with the ATcT value. [221 HO] At 
the W4 level, the calculated singlet-triplet splitting is found to be 609 cm -1 , somewhat lower 
than the experimental value [52] of 716 cm -1 (2.05 kcal/mol). Increasing the level of theory 
to W4.2 has a negligible effect on the triplet state (which has much weaker multireference 
character) but stabilizes the singlet state by a small but significant amount. As a result, at 
the W4.2 level the splitting goes up to 662 cm" 1 — the difference with W4 is entirely due to 
higher-order T 3 effects in the core-valence term. Further 'ramping up' to the W4.3 level sees 
a slight destabilization of the triplet state (mostly due to the T 3 — (T) and T 4 contributions 
being calculated with larger basis sets), but a somewhat more pronounced stabilization of 
the singlet state due to quintuple and sextuple excitations. Consequently, perfect agreement 
with experiment is reached, with an energy difference of 2.05 kcal/mol, or 2.00 kcal/mol 
if the experimental spin-orbit coupling constant [52], A=-15.25 cm -1 , is taken into account. 
(It is not entirely clear whether the experimental T e =716 cm -1 refers to the lowest spin- 
orbit component of the triplet state or to the average of the three spin-orbit components. 
Incidentally, our calculated A = —14.39 cm -1 using the same procedure as in Ref. [53].) 
For the isoelectronic BN diatomic — which exhibits even more pathological nondynamical 
correlation effects — we likewise found [53] that reproduction of the very small singlet-triplet 
splitting required accounting for connected sextuple and especially quintuple excitations, 
contributions of which have similar magnitudes as in C2. For the sake of completeness, 
our computed T e for the b 3 S~ state is 18.16 kcal/mol at the W4 level and 18.44 kcal/mol 
at the W4.3 level, compared to an experimental value [52] of 18.40 kcal/mol. Once again, 
essentially all post-W4 change is due to the pathological singlet state. 



Perusing both Tables |III| and |IV[ it becomes evident that higher-order correlation effects 
wax increasingly more prominent as n goes up. This phenomenon is illustrated graphically 
in Figure [l] For example, higher-order triples for C9 reduce TAE by 6.4 kcal/mol, while 
connected quadruples will increase it by approximately 9.8 kcal/mol. While these two effects 
partially compensate each other, together this still amounts to an n-particle truncation error 
at the CCSD(T) level of about 3.4 kcal/mol. For comparison, the same effect is just 0.7 
kcal/mol (W3.21ite) or 1.0 kcal/mol (W4) for C 3 , and 0.8 kcal/mol (W3.21ite) for C 6 . 

Connected quintuples likewise increase in importance as n grows, reaching nearly 1 
kJ/mol for C5 (the largest system for which we were able to calculate their contribution 
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explicitly). 

Interestingly enough, with the anomalous exception of C2( 1 £jj~), all molecules considered 
here have %TAE[(T)] diagnostics hovering in the 5% range, which indicates moderate but 
not severe nondynamical correlation |22j. 

We also note that the very small linear-cyclic energy difference in C4 has an appreciable 
post-CCSD(T) contribution: at the W4 level, higher-order connected triples and connected 
quadruples favor the linear structure by 0.5 and 0.4 kcal/mol, respectively. Likewise, in 
C6, higher-order triples favor the linear structure by 0.5 kcal/mol (W3.2) while connected 
quadruples only affect the equilibrium by 0.08 kcal/mol (favoring the ring). In polyacetylenic 
Cs, however, higher-order triples favor the ring by 2 kcal/mol, while the effect of connected 
quadruples is nearly an order of magnitude smaller. This raises the question whether a 
general trend exists: higher-order triples favoring polyacetylenic rings over linear structures 
while favoring linear over cumulenic ring structures. Unfortunately, CCSDT calculations 
beyond Cg are beyond our computational resources. 

The equilibrium structure of the C10 ring displays clear angle alternation. It was 
previously noted[6] that the barrier towards the -Dioh saddle point is very small, and there has 
been some speculation (e.g., Ref . [25 | l54"l 155] 156] ) that it might actually disappear at the basis 
set limit. Our calculated CCSD(T)/cc-pVTZ, CCSD(T)/cc-pVQZ(no g), and CCSD(T)/cc- 
pVQZ barriers are 1.11, 0.48, and 0.33 kcal/mol, respectively. At the {Q,5} basis set limit at 
CCSD(T)/cc-pVQZ optimum geometries, the D 5h alternating-angles ring is 23.00 kcal/mol 
more stable at the SCF level than the D W h saddle point. CCSD valence correlation lowers 
this by 18.37 kcal/mol to 4.63 kcal/mol, while parenthetical triples shave off yet another 
4.32 kcal/mol. These components add up to a basis set limit CCSD(T) deformation energy 
of 0.30 kcal/mol. MP2, as expected[6l [17], grossly favors the more symmetric structure, and 
SCS-MP2/cc-pVQZ[57j is unable to overcome this bias, reaching -26.75 kcal/mol for the 
isomerization energy. At the SCS-CCSD/ cc-pVQZ |58j level, however, +0.15 kcal/mol is 
obtained, slightly increasing to 0.18 kcal/mol at the SCS-CCSD/cc-pV5Z//CCSD(T)/cc- 
pVQZ level. These numbers are as close to the CCSD(T) answer as one can reasonably 
hope. 

Notwithstanding the above, the calculated B3LYP/cc-pVTZ zero-point vibrational 
energy of the saddle point is about 1 kcal/mol smaller than for the D 5 ^ minimum. This 
implies that the vibrationally averaged structure will have -D10/1 symmetry even at K. 
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Inner-shell correlation contributions become fairly hefty for the larger clusters, reaching 
12.6 kcal/mol for Cio- For the linear C n clusters, these contributions scale almost perfectly 
linearly with n. For Cq, Cg, and Cio, inner-shell correlation systematically favors the cyclic 
structures: the somewhat anomalous bicyclic C4 structure does not follow this trend. 

It was earlier noted by RB, for the linear clusters, that the reactions 2C n ^C n+ 2+C n _2 
are nearly thermoneutral at their level of theory. We do find the same to be the case at 
the W3.21ite level: however, relative to the (quite small) overall reaction energies, the post- 
CCSD(T) contributions are still non-negligible. If we assume that those reaction energies 
are converged at the W3.21ite level (which is probably justified to 0.5 kcal/mol, and certainly 
to 1 kcal/mol), we can obtain extrapolated W4 atomization energies for species heavier than 
C5 (the practical limit for full W4 calculations with the presently available hardware). Our 
best estimates are presented in Table [V) 

Perhaps the most striking feature of these numbers is the outstanding agreement between 
our best calculated TAE values for the linear clusters and the earlier, empirically scaled, 
data of MFG. Indeed, it is hard to believe that these latter rudimentary estimates - blissfully 
ignorant of higher-order correlation effects, inner-shell correlation, etc. — capture so much 
of the thermochemistry in these systems. We note that this scaling is not transferable to 
the rings, for which the MFG estimates are clearly too low. 

Our best values agree reasonably well with the available measured data of Gingerich and 
coworkers |2U]. considering the size of the uncertainties on the latter — although generally 
the calculated values are near the upper edges of the experimental error bars. 

Inspection of Table III also suggests (not surprisingly) that as n increases, the C n 
atomization energies at the W3.21ite level are progressively greater overestimates compared 
to the estimated W4 numbers. 

In conjunction with the ATcT revised heat of formation of carbon atom [23]^. 
Ai/j [C(g)] = 170.06±0.026 kcal/mol, we can offer the following revised heats of formation 
at K: C 2 ( 1 S+) 196.04, C 2 ( 3 n u ) 197.72, C 3 ( 1 S+) 194.33, C 4 ( 3 S;) 251.06, C 4 ( 1 A 9 ) 250.13, 
C 5 ( 1 S+) 253.63, C 6 ( 3 S 9 ) 304.14, C 6 ( 1 A , 1 ) 290.65, C 7 ( 1 S+) 312.94, C 8 ( 3 £-) 358.58, C 8 ( 1 A g ) 
345.47, C 9 ( 1 S+) 371.29, Ci ( 3 S") 412.33, and Ci ( 1 A' 1 ) 345.01 kcal/mol. 



Table VI presents computed atomization energies, and error statistics relative to our best 
values, for more approximate compound thermochemistry methods such as G2MP2[59J, 
G2jSD], G3jSI], G3B3[62], G4[53], G4(MP2)[M], CBS-QB3[ES], CBS-APNO[SB], WlhgU, 
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and W2.2 for the thermochemistry of these systems. 

Excluding C 8 and larger clusters (for which we deem our best calculated values to be less 
reliable), one sees similar error statistics for Wlh and for W2.2, suggesting that an RMSD of 
about 1 kcal/mol is the best one can hope for without including post-CCSD(T) correlation 
effects. CBS-APNO and Wlh have comparable RMSD errors: however, while the Wlh 
values tend to be underestimates (consistent with the positive post-CCSD(T) correlation 
effects), the CBS-APNO values tend to be overestimates. The next best performer, G4, 
errs on both sides. Somewhat counterintuitively, G2MP2 displays marginally better error 
statistics than G2 — despite MP2 being clearly inappropriate for these systems — and G3 
noticeably better ones than G3B3 — despite several qualitatively wrong reference geometries 
and spurious imaginary frequencies at the levels of theory used in G3 for reference geometries 
(MP2/6-31G*) and ZPVEs (HF/6-31G*), respectively. We believe the apparently better 
performance of G3 to be due to error compensation. CBS-QB3 clearly outperforms G3B3: 
in addition, it errs on the 'right' side — underestimating the best values, consistent with 
CCSD(T) limits. The unrealistically small linear-cyclic difference calculated for Cg is 
something of a 'clinker' on the part of CBS-QB3, however. One statement we can confidently 
make is that, for these systems too, G4 clearly outperforms the other members of the Gn 
series. Performance of G4(MP2), which actually does include a CCSD(T)/6-31G(d) step 
but considers basis set extension effects at the MP2 level only, is considerably degraded 
relative to G4, rather more so than generally seenjSl]. This is not surprising in light of the 
poor performance of MP2 for the present systems. 

Finally, we compare the relative performance of different DFT exchange-correlation 
functionals in predicting the binding energies of the carbon clusters considered in the present 
work. As reference data we use our best available |67j nonrelativistic, clamped-nuclei, zero- 
point exclusive TAEs. The exchange-correlation functionals employed include the following 
classes: (a) generalized gradient approximation (GGA): HCTH407,[68J BLYP, [Hll ED] 
BPW91,[69l[n] BP86[69l[72] andPBE[73j; (b) meta-GGA: M06-L,p3] VSXC,[75] TPSSpS] 
and r-HCTH;[77] (c) hybrid GGA: PBE0,[78j B97-2,[79j B3PW91,[7H EDj B97-1,(HI] 
B98,[82J TPSSh,[83] B3LYP, (TOj |80j [84] mPWlK|85l |86] and BHLYP;[87] (d) hybrid meta- 
GGA: mPWlB95,[86l|88] B1B95,[69J [89] M06,[90] PW6B95,[9T] TPSS1KCIS,[92] M05,[93] 
r-HCTHh,[77] BMK,|1] M06-2X,[11] BB1K[95J and PWB6K;[9TJ and (e) double hybrid 
functionals: B2-PLYP,[96J mPW2-PLYP,[97j B2T-PLYP, P] B2K-PLYPJSH] and B2GP- 
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PLYP.[99] Unless otherwise indicated, the pc-2 basis set of Jensen [TOO] , which is of 
[4s3p2dlf] quality but optimized for Hartree-Fock and DFT calculations, was employed 
throughout. The root mean square deviations (RMSD), mean signed deviations (MSD), 



and mean average deviations (MAD) are gathered in Table VII 

In general, the GGA, meta-GGA, hybrid-GGA, and hybrid meta GGA functionals 
systematically overestimate the binding energies and lead to unacceptable RMSDs of 11-57 
kcal/mol. While the so called 'kinetics' functionals (HF exchange > 40%) BHLYP, mPWIK, 
PWB6K, BB1K, and to a lesser extent BMK tend to systematically underestimate the 
binding energies. 

Not quite surprisingly, all conventional DFT functionals fare rather poorly, with RMSDs 
in the double digits overall. Simply eliminating the pathological C2 singlet case brings 
RMSD below 10 kcal/mol for four functionals: M05, BB1K, M06-2X, and BMK (9.4, 8.4, 
7.5, and 6.2 kcal/mol, respectively). Interestingly, once the linear systems larger than C4 
(and the C10 ring) are also eliminated, RMSDs drop sharply almost across the board. Long 
cumulenic chains are expected to be a 'torture test' for any DFT method. 

This becomes especially clear for the double hybrids. Near the basis set limit, RMSD 
for the reduced set of systems actually reaches the 2 kcal/mol range that can generally 
be expected for functionals like B2GP-PLYP[99]. For the long chains, however, severe 
overestimates set in. We note that these systems have increasingly narrower HOMO-LUMO 
gaps, and that [TUT] they have exceptionally large electron affinities for large n. It is perhaps 
not overly surprising that systems that put even post-CCSD(T) methods severely to the test 
(and, for instance, have connected quadruples contributions that climb to 10 kcal/mol for 
C9) would be beyond the reach of double hybrids, let alone conventional DFT functionals. 

In terms of singlet-triplet equilibria, as noted earlier. [T02] shortcomings of conventional 
DFT functionals (and particularly 'kinetics' functionals) are exacerbated by a natural bias 
of hybrid functionals towards high-spin states. This issue is mitigated by the MP2-like 
correlation in double hybrids. 

IV. CONCLUSIONS 

The thermochemistry of the carbon clusters C n (n=2-10) has been revisited by means of 
W4 theory and W3.21ite theory. Particularly the larger clusters exhibit very pronounced 
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post-CCSD(T) correlation effects. Despite this, our best calculated total atomization 
energies agree surprisingly well with 1991 estimates obtained from scaled CCD(ST)/6-31G* 
data. Accurately reproducing the small singlet-triplet splitting in C 2 requires inclusion 
of connected quintuple and sextuple excitations. Post-CCSD(T) correlation effects in C4 
stabilize the linear form. Linear/cyclic equilibria in C6, Cs, and C10 are not strongly 
affected by connected quadruples, but they are affected by higher-order triples, which favor 
polyacetylenic rings but disfavor cumulenic ones. Near the basis set limit, Ci does undergo 
bond angle alternation in the bottom-of-the-well structure, although it is expected to be 
absent in the vibrationally averaged structure. The thermochemistry of these systems, and 
particularly the longer linear chains, is a particularly difficult test for density functional 
methods. Particularly for the smaller chains and the rings, double-hybrid functionals 
clearly outperform convential DFT functionals for these systems. Among compound 
thermochemistry schemes, G4 clearly outperforms the other members of the Gn family. Our 
best estimates for total atomization energies at K are: C 2 ( 1 Sp ) 144.07, C 2 ( 3 II U ) 142.39, 
C 3 ( 1 S+) 315.83, C 4 ( 3 S 9 ) 429.16, C A { l A g ) 430.09, C 5 (^) 596.64, C 6 ( 3 £ 9 ) 717.19, C 6 ( 1 A / 1 ) 
729.68, C 7 ( 1 E+) 877.45, C 8 ( 3 £ 9 ) 1001.86, C^Ag) 1014.97, C 9 ( 1 E+) 1159.21, Ci ( 3 E 9 ) 
1288.22, and Ci (M / 1 ) 1355.54 kcal/mol. 
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TABLE I: Theoretical geometries (in Angstrom and degree). In the linear systems, bond lengths 
are numbered from the outside in. 
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1.296 


1.284 


1.271 






B3LYP/pc-2 


C 7 ( 1 E+ 




1.282 


1.284 


1.269 








c 8 ( 3 s- 

° V y 




1.291 


1.286 


1.270 


1.277 






Cq( 1 S + 

y \ y 


Dooh 


1.281 


1.286 


1.267 


1.273 






Cl()( 3 S g 


) Dooh 


1.288 


1.287 


1.269 


1.277 


1.270 




v y / 


D 2h 


1.442 








117.6 




CeC^i) 




1.319 








147.9 






D 4 fc 


1.252 


1.380 






162.1 






) D 5h 


1.288 








161.8 




C 2 ( 1 Sg h 


Doch 


1.2400 a 












c 2 ( 3 n„) 




1.3153 












C 2 ( 3 E" 




1.3728 












C3( 1 S+ 


d°° 


1.2981 












C4f 3 S~ 


^oon 


1.3140 


1.2936 










c 5 ( 1 sj 


Doch 


1.2936 


1.2857 










c 6 ( 3 s- 


Dooh 


1.3052 


1.2905 


1.2780 






CCSD(T)/cc-pVQZ 


c 7 ( 1 s+ 


Doch 


1.2933 


1.2901 


1.2759 








c 8 ( 3 s- 


Dooh 


1.3010 


1.2923 


1.2767 


1.2832 








Doch 


1.2930 


1.2928 


1.2745 


1.2799 






C 4 ( 1 A 9 ) 


D 2h 


1.4492 








117.09 






D 3 h 


1.3282 








148.79 






C4.f1 


1.2592 


1.3926 






162.74 




CioC 1 ^ 


) D 5h 


1.2940 








158.25 






) D 10 h 


1.2914 








[144.00] 



(a) Fixed reference geometry used for consistency with earlier post-W4 work [ID] . Actual CCSD(T)/cc-pVQZ bond distance is 
1.2458 A. W4 TAE e at that geometry is 143.87 kcal/mol, just 0.01 kcal/mol higher than at 1.24 A. 
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TABLE II: Diagnostics for importance of nondynamical correlation 





%TAE e %TAE e 


%TAE e 


%TAE e 






Largest T2 


NO occupations 




[SCF] a 


[(T)]» 


[post-CCSD(T)] a [T 4 + T 5 ] a 


diagnostic 


amplitudes 


HDOMO 6 


LUMO 
















- CCSD/cc-pVTZ - 




c 2 ( 1 s+) 


12.61 


13.29 


0.28 


1.79 


0.038 


0.086 


0.29 


1.629 


0.362 


c 2 ( 3 n u ) 


50.56 


6.51 


0.85 


1.02 


0.020 


0.039 


0.12 


1.934 


0.084 


c 2 ( 3 £-) 


67.70 


3.40 


0.53 


0.34 


0.011 


0.021 


0.08 (x2) 


1.933 


0.044 


c a ( 1 £+) 


64.88 


5.52 


0.30 


0.64 


0.023 


0.052 


0.10 (x2) 


1.913 (x2) 


0.074 (x2) 


c 4 ( 3 s-) 


66.09 


4.61 


0.30 


0.62 


0.018 


0.035 


0.07 (x2) 


1.915 (x2) 


0.068 (x2) 


C 4 ( 1 A 9 ) 


64.44 


4.70 


0.10 


0.54 


0.015 


0.034 


0.09 


1.908 


0.074 


c 5 ( 1 s+) 


66.64 


5.37 


0.14 


0.60 


0.024 


0.060 


0.09 (x2) 


1.903 (x2) 


0.084 (x2) 


C 6 ( 3 £-) 


66.73 


4.79 


0.10 


0.54 


0.023 


0.050 


0.08 


1.912 (x2) 


0.074 (x2) 




64.19 


5.00 


0.03 


0.52 


0.035 


0.116 


0.07 (x2) 


1.915 (x2) 


0.098 




67.18 


5.30 


0.19 


0.71 


0.025 


0.066 


0.09 (x2) 


1.897 (x2) 


0.090 (x2) 


c 8 ( 3 s-) 


67.22 


4.80 


0.20 


0.67 


0.025 


0.061 


0.08 


1.904 (x2) 


0.0866 (x2) 




66.15 


4.46 


0.43 


0.69 


0.039 


0.130 


0.06 (x2) 


1.889 d 


0.100 d 


C 9 ( 1 E+) 


67.59 


5.18 


0.29 


0.83 


0.025 


0.071 


0.08 (x2) 


1.895 (x2) 


0.096 (x2) 


Cio( 3 S-) 


67.58 


4.87 


N/A 


N/A 


0.026 


0.071 


0.05 (x2) 


1.901 (x2) 


0.907 (x2) 


CioC^iJDBfc 


67.95 


4.68 


N/A 


N/A 


0.037 


0.112 


0.07 (x2) 


1.900 (x2) c - d 


0.089 (x2) c ' d 



"Percentages of the total atomization energy relate to nonrelativistic, clamped-nuclei values with inner shell electrons 
constrained to be doubly occupied (C2-C5 from W4 theory, Cg and C7 from W3.2 theory, C8 and Cg from W3.21ite theory, 
C10 from W2.2 theory). 

b Highest doubly occupied molecular orbital 
c HOMO-1: 1.904(x2), LUMO+1: 0.075(x2) 
d cc-pVDZ basis set 
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TABLE III: Component breakdown of the final W2.2 and W3.21ite total atomization energies at 



the bottom of the well (in kcal/mol) obtained from B3LYP/cc-pVTZ reference geometries 





SCF 


valence valence T3 — (T) 


1 lxfO) 


inner 


rclativ. 


spin-orbit 


DBOC 


(a) 


TAE e 






CCSD 


CT1 
V 1 ) 




/cc-pVDZ 


shell 










W2.2 


W3 21ito 


C9f 1 E"h 


18.30 


1 07 fiO 


1 Q 

ly.oy 


-2 25 


2.95 


1.01 


-0.17 


-0.17 


0.03 


0.04 


lit: no 


1 A(\ 70 




73.36 


60.22 


9.19 


-0.17 


1.08 


0.96 


-0.09 


-0.17 


0.01 


0.03 


143.48 


144.41 


W V ^g ) 


86.79 


36.07 


4.27 


0.39 


0.32 


0.67 


-0.07 


-0.17 


0.00 


0.03 


127.56 


128.28 


^6 K ^g ) 


207.96 


93.01 


17.24 


-0.96 


1.65 


2.07 


-0.19 


-0.25 


0.02 


0.06 


319.87 


^90 




288.48 


125.46 


19.76 


-1.14 


1.99 


3.21 


-0.33 


-0.34 


-0.29 


0.08 


435.95 


436.84 


^ V y J 


281.76 


133.58 


20.08 


-1.74 


1.98 


2.41 


-0.29 


-0.34 


0.03 


0.08 


437.23 


437.51 


^0 V ^g J 


404.38 


167.41 


31.90 


-2.51 


4.01 


4.62 


-0.47 


-0.42 


0.07 


0.10 


607.50 


609.05 


*— 'D v 9 ) 


485.83 


204.97 


34.30 


-2.86 


4.20 


5.86 


-0.62 


-0.51 


-0.17 


0.12 


729.65 


731.06 


O V 1 / 


476.56 


227.68 


36.67 


-3.32 


4.13 


6.84 


-0.87 


-0.51 


0.06 


0.13 


746.43 


747.29 


C 7 ( 1 E+) 


599.93 


242.14 


46.43 


-4.04 


6.79 


7.22 


-0.76 


-0.59 


0.11 


0.19 


894.49 


897.34 


Cg( 3 S") 

o V ^9 ) 


683.35 


282.36 


48.82 


-4.79 


6.80 


8.50 


-0.92 


-0.68 


-0.10 


0.16 


1021.33 


1023.43 


Cs^Ag) 


681.28 


298.24 


45.97 


-2.62 


7.06 


9.53 


-1.16 


-0.68 


0.09 


0.17 


1033.28 


1037.80 


C 9 ( 1 S+) 


795.86 


317.23 


61.02 


-6.44 


9.79 


9.83 


-1.06 


-0.76 


0.15 


0.19 


1182.27 


1185.72 


Cl0( 3 Sg ) 


880.78 


359.06 


63.43 


N/A 


N/A 


11.13 


-1.21 


-0.85 


-0.02 


0.21 


1312.32 


N/A 


CioC 1 ^;) 


934.19 


376.34 


64.38 


N/A 


N/A 


12.60 


-1.46 


-0.85 


-0.01 


0.21 


1385.20 


N/A 




linear-cyclic isomerization energies 


C 4 ( 1 J 4„)^C 4 ( 3 S-) 


-6.72 


8.12 


0.32 


-0.60 


-0.01 


-0.80 


0.04 





0.32 


0.01 


1.29 


0.67 


Cef^'J-t-Ce^ET) 

"\ 1 / "\ y / 


-9.27 


22.71 


2.37 


-0.46 


-0.07 


0.98 


-0.25 





0.23 


0.00 


16.77 


16.24 


Cs^Ag^CsC^g) 


-2.07 


15.88 


-2.86 


2.17 


0.26 


1.03 


-0.24 





0.19 


0.00 


11.94 


14.37 


CioC^'J^CiofS-) 


53.41 


17.28 


0.95 


N/A 


N/A 


1.48 


-0.25 





0.01 


N/A 


72.88 


N/A 




isodesmic reaction energies 


2C5 ^C7+C3 


0.86 


-0.32 


0.12 


-0.03 


-0.42 


-0.05 


0.01 


0.00 


0.01 


-0.05 


0.62 


0.17 


2C6 ^C8+C4 


-0.17 


2.11 


0.02 


0.22 


-0.39 


0.01 


0.00 


0.00 


0.05 


0.00 


2.02 


1.85 


2C7 -^Cg+Cs 


-0.37 


-0.36 


-0.05 


0.87 


-0.22 


0.00 


0.00 


0.00 


0.01 


0.09 


-0.73 


-0.08 


2Cg — »Cl0+C6 


0.09 


0.70 


-0.09 


N/A 


N/A 


0.01 


0.00 


0.00 


-0.01 


0.00 


0.70 


N/A 


3C 5 ^C 9 +2C 3 


1.35 


-1.01 


0.20 


0.82 


-1.05 


-0.11 


0.02 


0.00 


0.02 


0.00 


0.48 


0.25 



(a) difference between the ACES II and MOLPRO definitions of the valence ROCCSD(T) 

(b) RCCSDTQ/cc-pVDZ-UCCSDTQ/cc-pVDZ 
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TABLE IV: Component breakdown of the final W4 C total atomization energies at the bottom of 



the well (in kcal/mol) obtained from CCSD(T)/cc-pVQZ reference geometries 







SCF 


valence valence T3 — (T) 


f 4 


i\ 


T 6 


inner 


rclativ. 


spin-orbit 


DBOC 


(a) 




TAE e a 








CCSD 


(T) 




(c) 






shell 














^2 V z-'g 1 




18.38 


1 07 fin 

1U / .UU 


1 Q ^7 


9 1 Q 


2.37 


0.24 


[01 


1.06 


-0.17 


-0.17 


0.03 


0.04 





146.52 


ditto W4.3 




18.38 


107.60 


19 37 


-2.24 


2.35 


0.32 


0.07 


1.25 


-0.17 


-0.17 


0.03 


0.00 





146.78 


CM 3 IL,~) 

^2 \ - L1 u j 




72.87 


60.65 


9.39 


-0.25 


1.33 


0.14 


[01 


0.89 


-0.09 


-0.17 


0.01 


0.03 


0.05 


144.78 


ditto W4.3 




72.87 


60.65 


9.39 


-0.27 


1.29 


0.12 


0.01 


0.92 


-0.09 


-0.17 


0.01 


0.00 


0.05 


144.72 


C9f 3 S~) 
^2 \ z-'g J 




86.62 


36.29 


4.35 


0.25 


0.41 


0.03 


roi 


0.64 


-0.07 


-0.17 


0.00 


0.03 


0.00 


128.36 


ditto W4.3 




86.62 


36.29 


4.35 


0.22 


0.42 


0.03 


0.001 


0.65 


-0.07 


-0.17 


0.00 


0.00 


0.00 


128.34 


Ca( 1 S+) 




207.18 


93.58 


17.63 


-1.10 


2.07 


-0.02 


[0] 


1.96 


-0.19 


-0.25 


0.02 


0.06 





320.90 


C 4 ( 3 T,-) 




287.62 


126.19 


20.06 


-1.37 


2.56 


0.12 


[0] 


3.08 


-0.34 


-0.34 


-0.35 


0.08 


0.33 


437.28 


C^Ag) 




280.97 


134.15 


20.48 


-1.91 


2.16 


0.19 


[0] 


2.28 


-0.29 


-0.34 


0.03 


0.08 





437.77 






403.09 


168.48 


32.47 


-2.75 


3.39 


Q.2\ b 


[0] 


4.42 


-0.48 


-0.42 


0.07 


0.10 





608.54 


c 6 ( 3 s-) c 




484.46 


206.00 


34.80 


-3.19 


3.92 


N/A 


[0] 


5.71 


-0.64 


-0.51 


-0.23 


0.12 


N/A 


715.89 


c 6 ( 1 ^i) c 




475.87 


228.16 


37.07 


-3.60 


3.85 


N/A 


[0] 


6.58 


-0.88 


-0.51 


0.06 


0.13 


N/A 


729.68 


C 7 ( 1 £+) c 




589.19 


243.37 


47.20 


-4.68 


6.37 


N/A 


[0] 


6.96 


-0.78 


-0.59 


0.12 


0.15 





896.22 




singlet-triplet and linear-cyclic energy differences 


G 2 (X 


^C 2 (a 3 n„) 


-54.49 


46.95 


9.98 


-1.94 


1.04 


0.10 


[0] 


0.17 


-0.08 





0.01 


0.01 


-0.05 


1.74 


ditto W4.3 




-54.49 


46.95 


9.98 


-1.96 


1.06 


0.20 


0.06 


0.33 


-0.08 





0.01 





-0.05 


2.05 d 


C 2 (X !£+)- 


-C 2 ( 3 S-) 


-68.24 


71.31 


15.02 


-2.44 


1.96 


0.21 


[0] 


0.41 


-0.10 





0.02 


0.01 


0.00 


18.16 


ditto W4.3 




-68.24 


71.31 


15.02 


-2.46 


1.93 


0.29 


0.07 


0.59 


-0.10 





0.02 





0.00 


18.44 e 


c 4 ( 1 a 9 )^c 4 ( 3 s-) 


-6.64 


7.96 


0.41 


-0.54 


-0.40 


0.07 


[0] 


-0.80 


0.05 





0.38 


0.01 


-0.33 


0.49 




-8.60 


22.16 


2.27 


-0.41 


-0.08 N/A 


[0] 


0.87 


0.24 





0.29 


0.00 N/A 


16.26 



(a) difference between the ACES II and MOLPRO definitions of the valence ROCCSD(T) 

(b) RCCSDTQ/cc-pVDZ-UCCSDTQ/cc-pVDZ 

(c) UHF reference: for RCCSDTQ-UCCSDTQ, cfr. (b) 
a Note that the TAE e do not include ADBOC 

b f 5 approx. as CCSDTQ(5) A /cc-pVDZ(no d)-CCSDTQ/cc-pVDZ(no d). 
C C B and C7 from W3.2 theory. 
^Experimental value: 2.05 kcal/moll52|. 
Experimental value: 18.40 kcal/mol l52l . 
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TABLE V: Total atomization energies at K (kcal/mol) a 





ZPVE 6 


W2.2 


W3.21ite 


W3.2 


W41ite 


W4 


MT 

cc-pVTZ cc-pVQZ 


MFG 


RB 


Exp. c 


uncert. 


c 2 ( 1 s+) 




2.64 d 


143.39 


144.06 


143.88 


143.93 


143.88 e 


145.1 


144.2 


144.8 


147.1 


144.6 


1.9 


c 2 ( 3 n„) 




2.34 d 


141.16 


142.07 


141.94 


141.99 


142.45-f 


141.1 


142.4 










c 2 ( 3 s-) 




2.09 d 


125.45 


126.19 


126.01 


126.13 


126.27 














c 3 ( 1 s+) 


lin 


5.07 


314.80 


315.52 


315.27 


315.32 


315.83 


314.4 


314.3 


315.7 


320.9 


311.4 


3.1 


c 4 ( 3 s-) 


lin 


8.11 


427.83 


428.72 


428.26 


428.34 


429.16 


423.5 


426.1 






428.5 


4.1 




eye 


7.68 


429.55 


429.83 


429.52 


429.60 


430.09 


427.7 


429.3 


430.3 


437.6 








lin 


11.89 


595.60 


597.16 


596.65 


596.81 


596.64 


592.9 


594.0 


596.8 


606.8 


591.5 


4.1 


C 6 ( 3 Ss) 


lin 


14.50 


715.15 


716.56 


715.89 


716. 01 9 


717.19 9 


708.0 


N/A 






706.3 


4.8 




eye 


16.98 


729.45 


730.31 


729.68 






720.6 


724.9 


722.8 


734.7 






C 7 ( 1 S+) 


lin 


18.18 


876.31 


879.16[878.8] 


878.03 


878.28' 1 


877.45 h 


870.9 


[873.7] 


877.7 


892.4 


873.1 


4.8 


c 8 ( 3 s s -) 


lin 


20.70 


1000.63 


1002.73[1004.4] 


1001.86* 






992.9 


N/A 






[984.0] 


[10.4] 


Cs{ x A g ) 


eye 


22.83 


1010.44 


1014.97 








998.9 


N/A 


1004.6 


1021.3 






C 9 ( 1 E+) 


lin 


24.35 


1157.91 


1161.37[1161.2] 


1160.36-? 


1160. 72^ 


1159. 2V 


1149.6 


N/A 


1158.4 


1177.8 


[1154.6] 


[10.4] 


Ci ( 3 £-) lin 


26.82 


1285.50 


1288.22'= [1288.9] 








[1277.8] 


N/A 










CioC 1 ^) 


eye 


29.65 


1355.54 










1343.1 


N/A 


1344.3 


1366.6 







Reaction energies (linear clusters) 



2C 4 -> 


c 6 +c a ( 3 n„) 


-0.80 


-1.18 


-1.32 




-2.10 






2C 4 -> 


C 6 +C 2 ( 3 £-) 


14.91 


14.70 


14.62 










2C 5 


^c 7 +c 3 


0.09 


-0.37 


0.01 




0.5 


-1.43 


9.9 


2C 6 


^C8+C4 


1.84 


1.67 






-0.4 






2C 7 


^c 9 +c 5 


-0.94 


-0.20 






-0.7 






2C 8 


^Cio+C6 


0.67 














3C 5 


^C 9 +2C 3 


-0.75 


-0.94 













"Reference geometries for all the W3.21ite calculations and for all clusters larger than C7 are at the B3LYP/cc-pVTZ level of 
theory; values in square parenthesis were obtained with the formula: TAEo[C n ]=2xTAEo[Cn— 2]— TAEo[C n — 4]. 
'The zero-point vibrational energies for C3 — C10 are B3LYP/pc-2 scaled by 0.985. 
c Ref. [20]. 

d From u>e/2 — u) e x e /i, with u) e and uj e x e taken from Ref. |52|. 

e ATcT value is 144.03±0.13 kcal/mol, W4.2 value is 144.05 kcal/mol (142.46 kcal/mol for the triplet state). 
W4.3 and W4.4 values for the singlet state are 144.08 and 144.07 kcal/mol, respectively HOI. 
f W4.3 value is 142.39 kcal/mol. 

9 Estimated by assuming that the 2C 4 ( 3 £g )^C 6 ( 3 Sg )+C 2 ( 3 II„) reaction energy remains unchanged at the W3.2 and 
post-W3.2 levels. Via C 2 ( 3 Sg ), the estimated W41ite and W4 numbers are 715.93 and 717.43 kcal/mol, respectively. 
^Estimated by assuming that the 2C5 — >C7+C3 reaction energy remains unchanged at the W3.2 and post-W3.2 levels. 
'Estimated by assuming that the 2C6 — >C8+C 4 reaction energy remains unchanged at the W3.21ite and W3.2 levels. 
•'Estimated by assuming that the 3C5 — >Cg+2C3 reaction energy remains unchanged at the W3.21ite and post-W3.21ite levels. 
fc Estimated by assuming that the 2C8 — >Cio+C6 reaction energy remains unchanged at the W2.2 and W3.21ite levels. 
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TABLE VI: Performance of standard composite computational thermochemistry methods. Total 



atomization energies at K and error statistics with respect to our best values (in kcal/mol) 





G2MP2 


G2 


G3 


G3B3 


CBS-QB3 CBS-APNO 


G4 


G4(MP2) 


Will 


W2.2 


Best estimatel67j 


c 2 ( 1 s+) 


146.81 


146.90 


146.94 


147.91 


144.29 


144.28 


147.01 


148.89 


143.15 


143.39 


144.07 


c 2 ( 3 n u ) 


140.28 


140.48 


142.72 


143.73 


140.13 


140.83 


142.35 


143.46 


141.00 


141.16 


142.39 


c 3 ( 1 s+) 


316.75 


316.89 


317.53 


318.86 


314.21 


314 96 


316.93 


318.84 


314.75 


314.79 


315.83 


C 4 (U g ) 


430.30 


430.73 


428.98 


430.91 


427.39 


430.12 


431.66 


433.94 


429.10 


429.52 


430.09 


C 4 ( a EJ) 


424.52 


425.63 


429.61 


431.65 


426.87 


428.10 


428.85 


429.67 


428.15 


427.73 


429.16 




596.95 


597.50 


599.98 


601.87 


594.52 


597.19 


598.34 


600.44 


595.71 


595.58 


596.64 




725.51 


726.01 


727.73 


732.43 


727.18 


731.10 


731.72 


733.92 


728.76 


729.36 


729.68 


C 6 ( 3 Sg ) 


710.69 


709.14 


718.39 


719.34 


714.86 


719 04 


715.19 


715.87 


715.44 


715.10 


717.19 


Ct( 1 S+) 


876.39 


877.42 


881.72 


884.02 


874.30 


878.83 


879.28 


881.41 


876.45 


876.28 


877.45 




1005.15 


1005.89 


1011.16 


1016.63 


1006.61 




1013.39 


1015.80 


1009.80 


1010.44 


1014.97 


c 8 ( 3 s-) 


994.13 


996.90 


1004.38 


1004.68 


1000.48 




999.09 


999.66 


1000.90 


1000.63 


1001.86 


c 9 ( 1 s+) 


1156.44 


1158.01 


1164.12 


1166.88 


1155.66 




1160.51 


1162.64 


1158.07 


1157.91 


1159.21 


CioC 1 ^) 


1350.14 


1351.63 


1359.15 


1363.49 


1353.80 




1359.14 


1361.25 


1355.11 


1355.54 


1355.54 


Cio( 3 Ss) 


1277.31 


1280.91 


1290.12 


1289.43 


1284.83 




1283.48 


1283.95 


1285.72 


1285.50 


1288.22 


RMSD 


5.39 


4.49 


2.79 


4.22 


3.22 




2.30 


3.45 


1.83 


1.76 




RMSD" 


4.65 


4.22 


2.77 


3.92 


3.30 




1.80 


3.11 


1.84 


1.73 




RMSD 6 


3.25 


3.41 


2.29 


3.58 


2.27 


1.16 


1.73 


3.31 


1.14 


1.17 




RMSD C 


3.31 


3.48 


2.21 


3.55 


2.41 


1.22 


1.51 


3.07 


1.17 


1.22 





a RMSD w/o Ci . 

6 RMSD w/o Cg-Cio. 

C RMSD w/o C 8 -Ci and C 2 ( 1 S+). 
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TABLE VII: Performance statistics (kcal/mol) of various exchange-correlation functionals for the 
carbon clusters considered in the present work. 0. Unless otherwise indicated all calculations were 



done with the pc-2 basis set. 



Class 


Functional 


RMSD 


MSD 


MAD 


RMSD° 


RMSD 6 




HCTH407 


32.7 


25.9 


27.8 


33.7 


15.3 




BLYP 


29.2 


22.0 


23.9 


30.1 


10.5 


GGA 


BPW91 


48.5 


40.5 


41.8 


50.2 


26.8 




BP86 


60.4 


51.6 


52.3 


62.7 


35.7 




PBE 


74.4 


63.9 


64.3 


77.2 


46.8 




M06-L 


43.3 


35.6 


37.4 


44.8 


23.5 


meta 


vsxc 


31.8 


24.8 


26.7 


32.8 


13.0 


GGA 


TPSS 


22.6 


15.6 


18.2 


22.9 


7.3 




tHCTH 


27.9 


21.2 


23.4 


28.6 


11.2 




PBEO 


28.1 


20.1 


23.9 


28.3 


14.1 




B97-2 


25.5 


18.6 


21.8 


25.7 


11.2 




B3PW91 


21.9 


14.3 


18.0 


21.6 


8.1 


hybrid 


B97-1 


22.1 


14.7 


18.3 


21.8 


7.8 


GGA 


B98 


14.5 


5.3 


11.7 


12.8 


5.5 




TPSSh 


12.2 


1.7 


9.7 


10.3 


7.1 




B3LYP 


13.1 


2.0 


10.9 


11.3 


8.2 




mPWIK 


20.2 


-17.7 


17.7 


17.0 


17.2 




BHLYP 


56.6 


-52.4 


52.4 


56.6 


48.1 




mPWlB95 


20.7 


13.6 


17.7 


20.3 


9.4 




B1B95 


18.1 


11.2 


15.1 


17.4 


7.3 




M06 


19.7 


14.2 


17.2 


19.6 


11.1 


hybrid 


PW6B95 


18.1 


10.9 


14.9 


17.4 


6.3 


meta 


TPSS1KCIS 


23.2 


16.2 


19.2 


23.3 


8.3 


GGA 


M05 


10.4 


4.3 


8.7 


9.4 


8.8 




i-HCTHh 


24.1 


16.8 


20.1 


24.2 


9,1 




BMK 


11.8 


-4.2 


7.5 


6.2 


7.6 




M06-2X 


10.5 


2.5 


7.7 


7.5 


1.2 




BB1K 


12.7 


-8.9 


9.0 


8.4 


9,9 




PWB6K 


14.6 


-11.1 


11.1 


10.4 


11.7 




B2GP-PLYP C 


9.1 


3.4 


7.0 


8.6 


2.6 




B2GP-PLYP d 


12.1 


7.1 


9.3 


12.1 


2.1 




B2GP-PLYP e 


7.7 


1.5 


6.0 


6.9 


3.8 




B2GP-PLYP/ 


10.2 


5.1 


7.8 


9.9 


1.8 




B2-PLYP C 


13.8 


8.7 


10.7 


14.0 


2.6 




B2-PLYP d 


16.4 


11.4 


12.9 


16.7 


4.1 




B2-PLYP e 


12.2 


7.1 


9.6 


12.2 


2.4 




B2-PLYP^ 


14.6 


9.7 


11.4 


14.8 


3.1 




B2T-PLYP C 


9.2 


3.3 


7.4 


8.6 


3.2 


double 


B2T-PLYP li 


11.7 


6.4 


9.1 


11.5 


2.0 


hybrid 


B2T-PLYP" 


8.0 


1.5 


6.5 


7.2 


4.3 




B2T-PLYP/ 


10.0 


4.6 


7.9 


9.6 


2.4 




B2K-PLYP C 


7.4 


0.4 


5.7 


6.4 


1.0 




B2K-PLYP d 


10.1 


4.8 


7.5 


9.8 


1.9 




B2K-PLYP e 


6.7 


-1.8 


5.2 


5.5 


5.5 




B2K-PLYP^ 


8.2 


2.5 


6.2 


7.5 


2.5 




mPW2-PLYP c 


8.1 


-0.1 


6.7 


6.9 


6.2 




mPW2-PLYP <i 


9.2 


2.4 


7.7 


8.3 


4.4 




mPW2-PLYP e 


7.8 


-1.6 


5.9 


6.3 


7.2 




mPW2-PLYP^ 


8.3 


0.8 


6.9 


7.2 


5.4 



"Excluding C 2 ( 1 S+). 

^Excluding C2( 1 'S^), Cio, and all the linear systems larger than C4. 
c (all electron) cc-pwCVQZ basis set. 

d (a\l electron) cc-pwCVQZ basis set combined with a CBS extrapolation where Nmin=15 as recommended in Ref. I103| . 
e (frozen core) pc-3 basis set. 

•f (frozen core) pc-3 basis set combined with a CBS extrapolation where Nmin=15 as recommended in Ref. |103| . 
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FIG. 1: T3— (T) and T4 contributions to the total atomization energies for the linear C n clusters 
(from W3.21ite theory, in kcal/mol) 
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